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ABSTRACT 

We report late-time spectroscopic observations of the Type Iln supernova (SN) 
1998S, taken 14 years after explosion using the Large Binocular Telescope. The optical 
spectrum exhibits strong broad emission features of [O i], [O n], and Ha, in addition 
to weaker features of [O m] , H/3 and [Fe n] . The last decade of evolution has exhibited 
a strengthening of the oxygen transitions relative to Ha, evidence that the late-time 
emission is powered by increasingly metal-rich SN ejecta crossing the reverse shock. 
The Ha luminosity of «8000 L Q requires that SN 1998S is still interacting with rela- 
tively dense circumstellar material (CSM), probably produced by the strong wind of 
a red supergiant progenitor at least ~ 10 3 years before explosion. The emission lines 
exhibit asymmetric blueshifted profiles, which implies that the receding hemisphere 
of the SN is obscured by dust. The [O in] A5007 line, in particular, exhibits a com- 
plete suppression of its red wing. This could be the result of the expected wavelength 
dependence for dust extinction or a smaller radial distribution for [O in]. In the lat- 
ter case, the red wing of [O in] could be absorbed by core dust, while both the blue 
and red wings are absorbed by dust within the cool dense shell between the forward 
and reverse shocks; this interpretation could explain why late-time [O in] emission 
from SNe is often weaker than models predict. The [O i] line exhibits double-peaked 
structure on top of the broader underlying profile, possibly due to emission from in- 
dividual clumps of ejecta or ring-like structures of metal-rich debris. The centroids of 
the peaks are blueshifted and lack a red counterpart. However, an archival spectrum 
obtained on day 1093 exhibits a third, redshifted peak, which we suspect has become 
extinguished by dust that formed over the last decade, after day 1093. This implies 
that the "missing" red components of multi-peaked oxygen profiles observed in other 
SNe might be obscured by varying degrees of dust extinction. 

Key words: supernovae: general — supernovae: individual (SN 1998S) 



1 INTRODUCTION 

Supernova (SN) 1998S, which occurred in the nearby host 
galaxy NGC 3877, is one of the most well-studied SNe of 
the core-collapse variety. It is designated a Type Iln, a 
class which has raised critical new questions about the fi- 
nal phases of massive-star evolution. The "Iln" designation 
refers to narrow emission lines of hydrogen observed in opti- 
cal spectra (see Fillipenko 1997), caused by relatively slow- 
moving circumstellar gas that is impacted by the blast wave 
or illuminated by UV radiation from the SN. The presence 
of these narrow lines requires dense circumstellar material 
(CSM) ejected by the progenitor star shortly (within a few 
years to decades) before the explosion. In some extreme 
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cases, up to 25 Mq has been ejected in the decade before 
core collapse (Smith et al. 2007, 2008, 2010; Woosley et al. 
2007), which has suggested a possible link to the class of 
eruptive luminous blue variables (LBVs). The impact of the 
SN blast wave on the CSM can efficiently convert kinetic 
energy into radiation, which explains why some Type Iln 
events have been associated with several of the most lumi- 
nous SNe yet recognized in the Universe, such as SN 2006gy 
(Smith et al. 2007, 2010; Ofek et al. 2007). In other less 
extreme cases, such as SN 1998S, SNe Iln can exhibit inter- 
action between the blast wave and a dense, but relatively- 
steady, hydrogen-rich wind from a red or blue supergiant 
progenitor. In all cases, long-term monitoring of the spec- 
tral evolution of Type Iln SNe present the opportunity to 
probe the geometry, kinematics, structure, and chemistry of 
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Figure 1. The local field of SN 1998S. The HST /WFPC2 image 
(F814W, left) is from 2000 Oct 1, and the LBT/MODS image 
(g-band; right) was obtained 2012 Jan 27. 



both the slow-moving CSM and the rapidly expanding SN 
ejecta. 

Early spectra of SN 1998S showed signs of CSM inter- 
action very soon after explosion. Typical of Type Iln evolu- 
tion at early times, relatively narrow Ha and He I emission 
lines exhibited broad Lorentzian wings, probably formed by 
electron scattering of photons diffusing out of the opaque 
pre-shock gas (Leonard et al. 2000; Chugai 2001; Fassia et 
al. 2001). Narrow P-Cygni profiles developed in the UV and 
optical spectra, generated by the slow-moving ionized stel- 
lar wind of the progenitor, which had an estimated velocity 
of 40-50 km s" 1 (Bowen et al. 2000; Fassia et al. 2001). At 
about 200 days, the Ha, H/3, and He I lines began to exhibit 
complex structure in the form of asymmetric multi-peaked 
profiles, interpreted as the result of an expanding disk-like 
structure within the outflow (Leonard et al. 2000; Gerardy 
et al. 2000). SN 1998S also showed signs of thermal dust 
emission in the IR, either produced by the SN itself, ther- 
mal reprocessing of SN radiation by CSM dust, light echoes, 
or some combination of these (Fassia et al. 2000; Gerardy 
et al. 2002). The SN was also a strong long-lasting source of 
thermal X-rays, and generated radio free-free emission after 
about 300 days (Pooley et al. 2002), additional signposts of 
a lively CSM interaction. 

An important question concerning SNe Iln is how 
closely the eruptive pre-SN mass loss is linked to the im- 
pending core-collapse itself, or for how long before explo- 
sion the progenitor exhibited an unusually high mass-loss 
rate. For Type Iln SNe, as well as for all other SNe that 
may explode into dense CSM, the impact of the blast wave 
upon the CSM can potentially manifest itself in observable 
characteristics that persist for decades. Such SNe, not nec- 
essarily classified specifically as Type Iln, include 1957D, 
1970G, 1979C, 1980K, 1986E, 1986J, 1993J, 1996cr, 1988Z, 
and of course SN 1987A (see Milisavljevic et al. 2012, and 
references therein). Here we present data obtained 14- years 
after explosion, which adds SN 1998S to this group of per- 
sistently interacting core-collapse SNe, and we discuss this 
object's emission-line features in context. 



2 OBSERVATIONS 

SN 1998S was targeted for spectroscopic observation on 2012 
Jan 27 (day 5079 after discovery) with the Multi-Object 
Double Spectrograph (MODS; Byard & O'Brien 2000) on 
the 8 m Large Binocular Telescope (LBT). Under ~0'.'8 see- 
ing, a source at the SN position was detected in the 60.0 sec 
r-band target acquisition image, shown in Figure [1] From 
this image we were able to extract a photometric measure- 
ment using aperture photometry, deriving an r-band (SDSS) 
brightness of 23.3±0.3 mag. 

Spectral images in the blue and red channels were ob- 
tained with MODS in longslit mode, utilizing the G670L 
grating and a l'.'2 slit, which yielded a spectral resolution of 
R «1000. Because of guide star constraints, we had to ob- 
serve the SN with the slit at ~45° from the parallactic angle, 
prompting the concern that some blue light could have been 
lost as a result. However, since the SN was observed at an 
airmass of only ~1.1, and a l'/2 slit was used under «0'.'8 
seeing, we do not expect substantial loss of blue flux. Four 
exposures totaling 80 min were obtained at the SN position. 

The spectral images were bias subtracted, flat fielded, 
and median combined with a suitable rejection filter to re- 
move cosmic rays. The spectrum was extracted using stan- 
dard IRAF routines. Removal of the background line emis- 
sion from the host galaxy required us to sample the back- 
ground from a narrow strip very close to the dispersion track 
of the SN emission, such that the background region over- 
lapped very slightly with the wings of the SN point spread 
function (PSF), but the result subtracted less than a few 
percent of the SN flux. The extracted spectrum was wave- 
length calibrated using spectra of HeNeAr lamps. The wave- 
length solution was corrected for the redshift of NGC 3877 
(z = 0.003, Tully & Fisher 1988). 



3 RESULTS 

The new LBT/MODS spectrum of SN 1998S from day 5079 
is shown in Figure [2] where it is compared to three earlier 
epochs from day 290, 474, and 1093 (archival spectra from 
Pozzo et al. 2004, obtained from the Weizmann interactive 
supernova data repositor^Q; Yaron & Gal- Yam 2012). For 
days 290 and 474, only the range of the Ha line is presented, 
since there was a lack of significant detections from other 
transitions at these epochs. Figure [3] shows our blue channel 
spectrum of SN 1998S and the earlier day 1093 epoch. For 
comparison, we included spectra of SN 1980K at 5509 days 
after discovery (Fesen et al. 1999; spectrum kindly provided 
by R. Fesen and D. Milisavljevic). Although the MODS 
wavelength coverage extends beyond the regions displayed 
in Figure 2, barely any continuum emission was detected 
on the latest epoch of SN 1998S outside of the displayed 
wavelength range. This insignificant amount of continuum 
is likely contaminated with background emission from the 
H n region that SN 1998S originated in, so we refrain from 
interpreting it. 

The red channel spectra are dominated by strong, broad 
emission lines of Ha, and the forbidden doublet transitions 
of [O i] AA6300, 6364 and [O n] AA7319,7330. There is also a 
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Figure 2. The LBT/MODS red channel spectrum of SN 1998S 
obtained 2012 Jan 27 (day 5079, black), compared with spectra of 
the same source on days 290, 474, and 1093 (red; Pozzo et al. 2004; 
Fassia et al. 2001). A spectrum of SN 1980K on day 5509 is also 
included for comparison (Fesen et al. 1999). On day 1093, the 
[O i] profile exhibited a triple-peaked narrow-component (black 
arrows), but the reddest of these has diminished by day 5079. 




4600 4800 5000 5200 5400 
Rest Wavelength (A) 

Figure 3. LBT/MODS blue channel spectrum of SN 1998S on 
day 5079 (black), compared to day 1093 (red, Pozzo et al. 2004), 
and SN 1980K 5509 days after discovery (orange, Fesen et al. 
1999) . 



al. 2001) and 15.5 ± 1.3 Mpc (NASA/IPAC Extragalactic 
Database [NED]), respectively. Extinction values were 
calculated for the rest wavelength of the emission line 
in question, assuming an extinction coefficient of i?=3.1 
(Cardelli, Clayton, & Mathis 1989). No correction was 
made for the uncertain degree of internal extinction. 

All of the line profiles exhibit asymmetric blueshifted 
peaks, probably the result of a suppression of flux on their 
red wings. The flux drop is particularly steep for [O in] and 
also for the central portion of the [O i] profile, with flux 
diminishing abruptly longward of their rest wavelengths. 



possible contribution from [Ca n] AA7291, 7324 which over- 
laps with [O n] , although the latter is likely to be the domi- 
nant transition (Fesen et al. 1999). A weak but broad feature 
of [Fe n] A7155 is also detected on the blue side of [O n], 
although [Ar ill] A7136 might be contributing to the flux as 
well. 

The blue channel is dominated by [O ill] A5007 and 
H/3 A4861, with the former exhibiting a dramatic cutoff of 
its blue wing. Weak [Fe n] emission is also present around 
A5 100-5200, where multiple transitions of this species lie. 

Emission-line intensities, luminosities, and velocity 
parameters are listed in Table 1. The full- width at zero 
intensity (FWZI) of the Ha line is ^240 A, or ^11,000 
km s _1 , while the oxygen lines have FWZI values closer to 
«8000-9000 km s" 1 . For H/3 and [O ill], the full extent of 
their respective red and blue wings is difficult to ascertain, 
because of the partial overlap of these lines. Luminosity 
calculations were made after applying foreground extinction 
and distance estimates of Ay = 0.68^0 25 ma g (Fassia et 



There are conspicuous changes between day 1093 and 
5079. Perhaps the most obvious are the disappearance of the 
intermediate- width, systemic- velocity components of the Ha 
line, and the increasing strength of oxygen emission with re- 
spect to Ha. Less obvious is the change in the triple-peaked 
profile of the [O 1] AA6300, 6364 emission feature. The cen- 
tral peak, assuming it traces the A6300 component of the 
doublet, appears blueshifted from rest by roughly —610 km 
s _1 , and the blue and red peaks are shifted by ~ —2600 
and +2800 km s _1 , respectively, from the central peak. It 
seems likely, although not certain, that each of these peaks 
represent different velocity components of the A6300 tran- 
sition of the doublet, since the underlying broad emission 
is centered near A6300, with no comparable broad emission 
centered near A6364. However, if the red peak is instead from 
the A6364 component of the doublet, then it is exhibiting a 
blueshift of —800 km s _1 . Although the change is subtle, the 
central peak appears to become more sharply cutoff on its 
red side by day 5079. The blueshifted peak remains more or 
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Table 1. Lino Measurements for SN 1998S on day 5079 (ssl4 yr). Luminosity errors are ssO.15 dex, mainly dependent on distance and 
extinction uncertainties (day 2148 error is 0.25 dex, see text). Line ratios for the spectroscopically similar SNe 1980K and 1979C, at ages 
of 14—15 yr and 12—14 yr, respectively, are included for comparison (Fesen et al. 1999). Approximate luminosity ratios are given with 
respect to Ha. CF94 model predictions are included in the last column for power-law, and lower limits for RSG models. 
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Figure 4. Ha luminosity curves for SN 1998S and SN 1980K. The 
luminosity at epochs <1000 days is consistent with the 111.5-day 
mean lifetime of radioactive 56 Co (dashed line), but later epochs 
from day 1093 and onward fall above this trend by as much as 10 
dex, requiring CSM interaction. The late-time luminosity of SN 
1998S appears most consistent with model predictions for an RSG 
density profile (CF94; solid line), while SN 1980K lies between 
this and a power-law density profile (dotted line), as discussed in 



less unchanged. Remarkably, the reddest peak is no longer 
detected on day 5079. 

The most recent epoch on day 5079 exhibits 
L(Hq) = (3.1 ± 1.4) x 10 37 erg s" 1 (f»8000 L©). To compare 
this with past luminosity values, we adopted flux measure- 
ments from Fassia et al. (2001) and Pozzo et al. (2004) from 
day 1093 and earlier, and converted them to luminosities 
using the same distance and extinction values as before. 
In addition, an Ha flux for day 2148 was extracted from 



a spectrum presented in Fransson et al. (2005; kindly pro- 
vided by T. Matheson, but not shown here). The day 2148 
spectrum was observed under poor seeing conditions. Thus, 
there was substantial and inseparable contamination from 
the background H n region, specifically, narrow Ha and [N 
n]-doublet emission. We interpolated over these narrow fea- 
tures in order to estimate the broad underlying Ha flux of 
the SN. As a result, we may have removed some fraction of 
narrow Ha flux potentially intrinsic to the SN, if any were 
present at rest velocity on day 2148. In light of these caveats, 
we estimate a conservative error of 0.25 dex for the luminos- 
ity on day 2148. For additional comparison, we also include 
the Ha light curve for SN 1980K in Figure g] (See Milisavl- 
jevic et al. 2012, and references therein), since this SN has 
exhibited spectral evolution very similar to SN 1998S, and 
has excellent late-time coverage. For SN 1980K, we adopted 
E(B — V) = 0.41 mag (Fesen et al. 1999) and a distance of 
6 Mpc (NED). 

Earlier than 1000 days, the light curve for Ha is well 
matched by the radioactive decay trend of 56 Co (111.5- 
day mean lifetime; Firestone et al. 1996). After 1000 days, 
the light curve departs from this trend, and hits a floor of 
~ lO 37 ' 5 I/0 that is maintained through day 5079. Since the 
day 5079 spectrum differs from the spectrum at early times, 
the late time Ha luminosity is not likely to be from an echo. 
Rather, the luminosity is accounted for by SN/CSM inter- 
action, which we discuss in the following section. 



4 DISCUSSION 

4.1 Source of the Broad Line Emission 

The collision between the SN explosion and the CSM can 
have lasting effects on the optical emission-line spectrum, 
long after the continuum has faded away. The impact of the 
SN blast wave on the CSM leads to the formation of a for- 
ward shock that meets the progenitor wind, and a reverse 
shock that gets bombarded by the expanding metal-rich SN 
ejecta (Chugai 1991; Chevalier & Fransson 1994, hereafter 
CF94). Gas in the cool dense shell (CDS) that develops be- 
tween the forward and reverse shocks absorbs far-UV and 
X-ray photons generated within the shock, and this excites 
low-ionization transitions such as [Fe n] and [O i] . Inside the 
reverse shock, the hot outer layers of the metal-rich ejecta 
are also ionized, but give rise to higher ionization transi- 
tions, such as [O ill] and UV metal lines. Both the CDS and 
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the ionized ejecta regions produce comparable amounts of 
Ha, H/3, and [O i] emission, while [O ill] is generated almost 
exclusively on the skin of the ejecta, according to CF94. At 
late times (t > 10 yr), the luminosity of the reverse shock 
becomes significantly higher than the forward-shocked CSM 
gas (by a factor of ~30, CF94). As the ejecta expands, [O i], 
[O n], and [O in] emission are expected to increase relative 
to Ha, as collisional de-excitation becomes less effective for 
these forbidden transitions as the density drops. 

CF94 provide predictions for the UV and optical 
emission- line luminosity that results from SN/CSM inter- 
action. The luminosity of a given line depends on the ioniz- 
ing luminosity of the reverse shock, which is sensitive to the 
density profile of the expanding supernova gas. This profile 
is, in turn, dependent on the stellar structure, the wind of 
the progenitor, and on the shock kinematics. CF94 consider 
the density profile of a red supergiant (RSG), which has rel- 
atively steep density gradients, in addition to more smooth 
power-law forms. The luminosity dependence of Ha can thus 
be expressed as 

L(Ha) oc — oc — — , (1) 

is Is \V W J 

where L 3 , T s , and V s are the shock ionizing luminosity, 
temperature, and velocity, respectively. M/v w is the wind 
density parameter, where the mass-loss rate and progenitor 
wind velocity are given in units of 1O -5 M0 yr -1 and 10 km 
s _1 , respectively. 

The Ha light curve for SN 1998S is shown in Figure 
3] and includes model tracks for the SN density profile of a 
RSG and power-law distribution. The power law curve con- 
nects epochs 2, 5, 10, and 17 yr, adopted directly from CF94 
(their Table 6); and the RSG curve was generated for epochs 
3, 5, 10, and 14 yr by scaling the 10-yr Ha luminosity given 
in CF94 by the time- variable L s /T s ratio. Both models have 
M/v-w — 5, which is an appropriate wind density parameter 
for SN 1998S, since the mass- loss rate and wind velocity of 
the progenitor have been estimated at 40-50 km s" 1 and 
l-2xl0" 4 M Q yr -1 (Fassia et al. 2001; Pooley et al. 2002). 
The RSG and power-law CF94 models have ejecta veloci- 
ties of 4318 and 5038 km s , respectively. These velocities 
imply that 5079 days after explosion the shock should be lo- 
cated ~2x 10 12 km from the SN origin (~13000 AU). Red 
supergiant winds extend to at least ~ 3 x 10 13 km radius 
(Lloyd, O'Brien & Kahn 1995). Thus, the ongoing SN/CSM 
interaction during the most recent epoch is to be expected. 

Overall, the late-time Ha luminosity trend of SN 1998S 
illustrated in Figure U is apparently most consistent with 
the RSG model. The Ha luminosity of SN 1980K, on the 
other hand, has late-time luminosity values «6 times lower 
than SN 1998S at similar epochs. The lower values could be 
explained as the result of the lower mass-loss rate of the SN 
1980K progenitor (2.5 x 10~ 5 Af© yr" 1 , Fesen et al. 1999), 
which is a factor of ~6 lower than SN 1998S. As a result, 
the wind density parameter and, hence, the reverse shock 
luminosity for this object would also be lower by a com- 
parable factor. Alternatively, the power-law density profile 
may be more consistent with SN 1980K, as argued by Fesen 
et al. (1999). However, although the RSG model is favored in 
the case of SN 1998S, it would premature to overly discrim- 
inate between either models using line luminosities alone, 
since one only has to tweak the wind density parameter or 



shock velocity to obtain a reasonable match to the data. It 
is also important to keep in mind that the luminosity could 
be affected by the highly uncertain degree of local extinc- 
tion and/or asymmetries in emitting volume. If the ejecta or 
progenitor wind are asymmetric, the models we have consid- 
ered here would be providing an underestimate of M/v w and 
L(Ha). In any case, the implied mass loss rate of ~ 10 _4 Mq 
yr" 1 for SN 1998S is extremely high, at the far upper end 
of the range exhibited by RSGs (Gehrz & Woolf 1971; Jura 
& Kleinmann 1990). 



Elevated mass-loss rates could be a sign of impending 
core collapse for RSGs, perhaps driven by pulsation insta- 
bilities that generate a superwind (Yoon & Cantiello 2010). 
A suitable Galactic analog to the progenitor of SN 1998S 
could be the extraordinary object VY CMa (see Smith et 
al. 2001 and references therein). This RSG star is experi- 
encing a very large mass-loss rate of M = 2-4 x 10~ 4 (with 
possible brief bouts as high as ~ 10" 3 Mq yr" 1 ). It exhibits 
an average wind speed of 30-40 km s~ and is surrounded 
by dense CSM that extends ~1000 AU from the star, esti- 
mated to have been ejected over the last ~1000 yrs. For these 
reasons, VY CMa has been highlighted as a viable SN Iln 
progenitor (Smith, Hinkle, & Ryde 2009a). For the progeni- 
tor of SN 1998S, the wind speed of 40-50 km s" 1 combined 
with our estimated shock distance of 13000 AU implies that 
its high rate of mass loss began at least ~1200 years before 
explosion. If this is close to the onset time for RSG super- 
winds (as suggested by Yoon & Cantiello 2010), and if the 
entering of this phase signifies a countdown to core collapse, 
it thus implies that VY CMa (whose superwind has been 
blowing for ~1000 years) could explode relatively soon. 



We observe an increase in oxygen emission relative to 
Ha over the last decade, providing evidence that the late- 
time emission traces increasingly metal-rich SN ejecta cross- 
ing the reverse shock (CF94). Between days 1093 and 5079 
the [O l]/Ha ratio has increased from 0.1 to 0.6. This is 
roughly consistent with the power-law model and, to a lesser 
extent, the RSG model. However, over the same time period 
the [O ill] /Ha ratio increased by a much smaller amount, 
from 0.2 to 0.3, which is an order of magnitude below the 
value predicted for the power-law model, and a factor of 
at least ~17 below the RSG model. Serious discrepancies 
between the observed late-time line ratios and the CF94 
model predictions, particularly regarding [O ill], have been 
reported for several other SNe observed at late times, in- 
cluding SN 1980K and 1979C (CF94; Fesen et al. 1999). 
The ratio is highly sensitive to the exact value of the shock 
velocity and the ejecta density. CF94 have suggested that a 
flatter density profile for the high- velocity ejecta, in addition 
to variations in shock temperature, could result in a lower 
ratio that would help resolve the discrepancy. Alternatively, 
the structure of the [O in] line profiles in SN 1998S and 
1980K could be an indication that the optical depth of the 
[O ill] emitting region beneath the reverse shock is larger 
than expected, which would also result in a lower observed 
line ratio; and we explore this possibility further below. The 
geometry and filling factor of the [O ill] emitting region(s) 
could affect the result as well. 
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Figure 5. Velocity profiles for the strongest lines in the spectrum 
of SN 1998S. The shortest wavelength transitions, particularly [O 
i] A6300 and [O in] A5007, exhibit the most extreme flux cutoff on 
the red side of the line at zero velocity, strong evidence that inter- 
nal dust extinction is the cause of the blue-shifted profiles. Note, 
the velocity of [0 n] + [Ca n] should be taken with caution, since 
the center of the line blend could not be accurately determined. 



4.2 Line Profiles and Dust Extinction 

The emission lines in the spectra SN 1998S exhibit sub- 
stantial asymmetry and a net blueshift of their centroids. A 
commonly proposed explanation for this phenomenon is ex- 
tinction by dust. In the case of SNe that interact with dense 
CSM, dust is expected to form in the cool dense shell (CDS) 
that develops between the forward and reverse shocks. Emis- 
sion lines formed in the post-shock gas, exterior to the dusty 
CDS, will suffer little or no extinction, while the same emis- 
sion features emerging from the receding hemisphere will be 
obscured by the dust. This results in the suppression of the 
red wings of the emission lines, causing a net blueshift. Dust 
formation in the CDS has been postulated to explain the 
spectral evolution of the Type II SNe 1980K, 1979C, 1993J, 
and 1999em (Dwek et al. 1983; Fesen et al. 1999; Wang & Hu 
1994; Elmhamdi et al. 2003), the Type Iln SNe 2005ip and 
2010jl (Smith et al. 2009b, 2012), and the Type Ib/Ibn SN 
2006jc (Smith et al. 2008). Blueshifted emission lines have 
also been observed to develop in the spectra of core-collapse 
SNe that do not exhibit evidence for immediate interaction 
with dense CSM. In fact, SN 1987A was the first case to 
provide a convincing connection between the development of 
emission-line asymmetries and dust formation. After about 
day 500, a drop in optical flux was accompanied by a rise 
in IR emission; while around the same time, emission lines 
from SN 1987A began to exhibit asymmetry and blueshift 



Figure 6. Ratio of the red-wing to blue-wing flux plotted against 
wavelength for SN 1998S (circles) and SN 1980K (triangles). The 
[O III] value for SN 1998S could be taken as an upper limit, since 
its red wing appears indistinguishable from pure continuum. The 
inset illustrates the blue and red wing regions of the Ha line of 
SN 1998S, centered on the rest wavelength. 

(Danziger et al. 1989; Lucy et al. 1989; Colgan et al. 1994; 
Wang et al. 1996). Similar phenomena was observed in the 
photometric and spectral evolution of SN 19901 (Elmhamdi 
et al. 2004). In these cases, however, dust formation likely 
occurred within the metal rich ejecta. 

A large dust component has already been postulated in 
the case of SN 1998S. Gerardy et al. (2000) and Pozzo et al. 
(2004) monitored the NIR photometry and spectra between 
days 305-1242, observing the onset of CO emission and the 
development of a significant IR excess in the spectral energy 
distribution. These features could be attributable to thermal 
emission from dust synthesized in the SN, or IR echoes re- 
flecting off of pre-existing dust clouds in the outer CSM (see 
Sugerman et al. 2012), or both. However, evolving asymme- 
try and increasing blueshift of the Ha and He I A10830 lines 
observed by Pozzo et al. (2004) provided evidence that dust 
is likely being synthesized within the SN CDS. 

The emission-line profiles of SN 1998S provide strong 
evidence for internal dust. Figure [5] provides a close-up of 
the most prominent lines plotted in velocity coordinates. 
All of the emission lines exhibit significant blueshift and ap- 
pear suppressed on their red wings. Interestingly, a relatively 
abrupt drop in flux at zero velocity is observed for the two 
shortest wavelength transitions of oxygen. The effect is most 
extreme for [O ill] A5007, which exhibits an apparently com- 
plete suppression of its red wing. Figure O shows the ratio 
of red-wing to blue- wing flux for [O ill], [O i], and Ha (H/3 
was excluded since [O ill] overlaps with its red wing, and 
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[O ll] + [Ca n] was excluded because the central wavelength 
of the uncertain line blend could not be accurately deter- 
mined). The trend could possibly be the result of wavelength 
dependent extinction. The flux drop between Ha and [O ill] 
could imply a relative extinction ratio of ~1.4 mag between 
A5007 and A6563, which is consistent with an 7? « 3 extinc- 
tion coefficient for the obscuring material (Cardelli, Clayton, 
& Mathis 1989). However, we stress caution in drawing any 
definitive quantitative conclusions about the extinction law 
from the line profiles. In fact, the relative blue- and red- 
wing flux ratios from one atomic transition to another could 
be affected by differences in the radial distributions of the 
[O i], [O in], and Ha emission zones. We discuss this further 
below. 



4-2.1 Evidence for Dust in the Unshocked Ejecta 

The relative weakness of [O ill] emission, and the complete 
suppression of its red wing imply some interesting possibili- 
ties for the distribution of obscuring dust within SN 1998S. 
Since the [O ill] emission is expected to be generated ex- 
clusively on the outer skin of the unshocked ejecta (CF94), 
interior to the CDS and reverse shock, then we expect dust 
formed in the CDS will obscure both sides of this emission 
line equally. Thus, the large net blueshift of the [O ill] emis- 
sion line and the lack of a red wing require the obscuring 
material to be located interior to the [O III] zone, within 
the unshocked SN ejecta. If so, then practically the entire 
receding hemisphere of the [O ill] emission skin would be 
obscured by dust in both the inner ejecta and in the CDS, 
explaining the lack of red wing flux. By comparison, the re- 
ceding components of Ha and [O i] emission, both of which 
can be generated at larger radii in the post-shock gas, and 
also potentially echoed by the CSM, would suffer less com- 
plete obscuration from inner ejecta dust, resulting in more 
red-wing flux for these lines relative to [O ill]. Since both the 
blue and red wings of [O ill] are absorbed by dust within the 
CDS, this could explain why late-time [O ill] emission from 
SNe such as SN 1998S, 1980K, and 1979C, is often weaker 
than models predict (CF94; Fesen et al. 1999). 

The spectrum of SN 1998S thus provides compelling 
evidence for dust formation in the CDS and within the 
unshocked ejecta. Dust formation in the latter zone has 
also been postulated in the cases of SNe 2005ip and 2006jc 
(Smith et al. 2008, 2009b). In addition, mid-infrared obser- 
vations of the SN remnant IE 0102.2-7219 in the Small 
Magellanic Cloud (Rho et al. 2009; Sandstrom et al. 2009) 
have provided evidence for the formation of dust within un- 
shocked SN ejecta, interior to the [O ill] emission zones. 

4.3 Evolution of Multi-peaked Oxygen Profiles 

The oxygen emission lines in the spectrum of SN 1998S ex- 
hibit interesting higher-order structure. Figure shows a 
close-up of the [O i] emission feature on days 1093 and 5079. 
On day 1093, the line exhibits three narrow peaks with max- 
ima at A6233, A6287, and A6346, all on top of a broader un- 
derlying profile. By day 5079, the reddest of these peaks is no 
longer detected, leaving behind a double-peaked asymmetric 
profile having an overall blueshift, similar in appearance to 
SN 1980K on day 5509. As noted in §3, although we cannot 
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Figure 7. Close-up of the [O i] A6300, 6364 profile on days 1093 
and 5079. The dashed line marks the rest wavelength of the A6300 
component of the transition. The dotted lines mark the center 
wavelengths/velocities of the triple peaks, assuming they all rep- 
resent A6300 (see caveats in text). SN 1980K on day 5509 exhibit 
very similar peak structure. 

completely rule out the possibility that the reddest peak on 
day 1093 is a blueshifted A6364 component of the doublet, 
we suspect that it actually represents a redshifted A6300 
component. Assuming this to be the case, the blue and red 
peaks would have velocities of approximately ±2600-2700 
km s _1 with respect to the central peak. The origin of these 
peaks and the cause of their temporal change is puzzling. 

Double-peaked oxygen emission features have been ob- 
served in many other core-collapse SNe (e.g., see Modjaz et 
al. 2008, and references therein; Milisavljevic et al. 2009). 
One interpretation involves emission from large individual 
clumps of oxygen-rich ejecta (Wang & Chevalier 2002). How- 
ever, this model could have difficulty explaining the appar- 
ent similarity of the wavelength spacing between peaks that 
is typically seen from one SN to another. Other interpreta- 
tions include emission from the equatorial ring or torus of a 
bipolar outflow viewed nearly edge-on (Chugai & Danziger 
1994; Maeda et al. 2008). However, the toroidal model al- 
ways produces two peaks centered around the rest velocity, 
in addition to a possible central component. Milisavljevic et 
al. (2010) pointed out that the double-peaked oxygen pro- 
files observed in core-collapse SNe (mainly Type lb, lib, and 
Ic) typically come in two main groups: symmetric profiles, 
where the trough between two peaks lies near rest velocity, 
with the peak separation lying conspicuously close the the 
64 A separation of the doublet components; and asymmetric 
profiles, where one component lies near rest and the other is 
blueshifted, having no conspicuous redshifted counterpart. 
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The [O i] profile of SN 1998S on day 5079 falls into the 
latter asymmetric category. On day 1093, however, the line 
appears more symmetric with the inclusion of the conspic- 
uous red peak. The subsequent disappearance of this peak 
by day 5079 implies that it has been extinguished by dust, 
and therefore, originates from the receding hemisphere of 
the SN. This supports the interpretation that the peak is a 
redshifted A6300 emission feature from the doublet, and not 
a blueshifted A6364 feature. 

We question whether the [O i] profile evolution could 
be related to the emergence and disappearance of the triple- 
peaked Ha profile at earlier times. By day 258, the Ha profile 
of SN 1998S exhibited a dramatic triple-peaked profile (see 
Figure [2]). Its red and blue peaks were shifted by roughly 
±4000 km s _1 with respect to the central peak. Spectra 
from subsequent epochs exhibited a continuous flux decrease 
of the red peak, attributed to increasing extinction of the 
receding hemisphere of the SN by dust. By day «650 the 
red peak had become extinguished below detectability. The 
evolving red and blue peaks have been interpreted as the 
result of emission from an expanding ring or torus, with 
the central peak possibly arising from shocked circumstellar 
clumps (Gerardy et al. 2000; Leonard et al. 2000). How- 
ever, Fransson et al. (2005) identified inconsistencies be- 
tween this model and the observations, and instead favored 
a scenario in which double-peaked emission profiles emerge 
from a geometrically thin CDS behind the reverse shock, 
at large optical depth. Still, the CDS would probably have 
to be non-spherical to be consistent with the polarization 
measurements by Leonard et al. (2000). In any case, the ob- 
served evolution of the triple-peaked Ha profile motivates 
the following question. Could the triple-peaked structure of 
[O i] A6300 seen on day 1093 — and the subsequent dis- 
appearance of its red component by day 5079 — be at- 
tributable to similar processes involving ring- like/toroidal 
ejecta or a non-spherical thin CDS, where the red peak of a 
corresponding triple-peaked profile eventually becomes ex- 
tinguished as a result of persistent dust formation? Indeed, 
whatever structure the reverse-shock/CSM interaction ex- 
hibited at relatively early times would likely manifest again 
later as oxygen-rich ejecta eventually crosses over. 

Disk, ring, and toroidal structures in SN ejecta have 
been suggested as viable explanations for observations of 
oxygen-rich SN remnants such as SNR IE 0102.2-7219 in 
the Small Magellanic Cloud and Cassiopeia A in the Galaxy, 
both of which exhibit velocity structures indicative of asym- 
metric explosion geometries (Tuoho & Dopita 1983); indeed, 
the integrated spectra of the [O i] A6300 line in Cassiopeia 
A exhibits a multi-peaked structure similar to what is seen 
in the spectra SN 1998S and 1980K (Milisavljevic et al. 
2012). In addition, spectropolarimetry of other core-collapse 
SNe (Wang & Wheeler 2008, and references therein), the 
morphologies of many young Galactic remnants, and pulsar 
space velocity distributions (Cordes, Romani, & Lundgren 
1993; Lyne & Lorimer 1994; Cordes & Chernoff 1998), col- 
lectively support the commonality of asymmetric explosions 
from core-collapse SNe. 



5 CONCLUSIONS 

Our observations of the Iln SN 1998S show that the CSM 
interaction is still going strong 14 years after explosion. 
The Ha luminosity has hit a slowly declining floor of L = 
3.1 x 10 37 erg s' 1 (^8000 L Q ). This behavior is very sim- 
ilar to SN 1980K, but indicates that the pre-SN mass-loss 
rate for the RSG progenitor of SN 1998S was extraordinarily 
high, even «1200 years before core collapse. This behavior 
suggests that other extreme RSGs, such as VY CMa, can 
potentially endure a similarly high degree of mass loss for 
at least 1200 years before exploding. 

The blueshifted emission-line profiles observed from SN 
1998S provide compelling evidence for the formation of dust. 
The relatively strong extinction of the red wing of the [O ill] 
A5007 line, in particular, is very interesting, as this is either 
a consequence of the wavelength dependence of the extinc- 
tion, or the fact that the [O in] emission is generated exclu- 
sively interior to the shocks on the ionized skin of the inner 
ejecta, having its receding hemisphere completely obscured 
by inner core dust. Evidence for dust within the unshocked 
ejecta of SNe, interior to the [O in] emission zone, was ac- 
quired from mid-infrared observations of the SN remnant 
IE 0102.2-7219 in the Small Magellanic Cloud (Rho et al. 
2009; Sandstrom et al. 2009). Dusty O-rich remnants such as 
this could thus provide a glimpse at the possible appearance 
of SN 1998S and related SNe, centuries into the future. 

Finally, the observed evolution of the multi-peaked [O i] 
profile, particularly the disappearance of the individual red 
peak, suggests that the red peak originated within the re- 
ceding hemisphere of the SN and became obscured by dust 
that formed between days 1093 and 5079. The peak would 
therefore have to represent the A6300 component of the dou- 
blet. This gives some credence to the interpretation that the 
individual peaks of oxygen profiles observed in other SNe 
represent emission features from separate structures within 
the outflow. This also implies that the "missing" conspicu- 
ous red peak in the multi-peaked [O i] profiles from other 
SNe might be the result of varying degrees of dust obscura- 
tion. Interestingly, SN 1998S must have formed significant 
amounts of dust after 3 years in order to extinguish the red 
peak between days 1093 and 5079. 
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